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Abstract 



. A solution to the cosmological constant problem has been proposed in which 

CN ■ our universe is a 3-brane in a 5-dimensional spacetime. With a bulk scalar, 

the field equations admit a Poincare invariant brane solution regardless of 

O - 

. the value of the cosmological constant (tension) on the brane. However, the 

o , 

solution does not include matter in the brane. We find new exact static 

O 

solutions with matter density and pressure in the brane. We study small 
Qh. perturbations about these solutions. None seem consistent with observational 



cosmology. As a byproduct we find a class of new matterless static solutions 



^ . and a non-static solution which, curiously, requires the string value for the 

SH ; 

dilaton coupling. 



I. INTRODUCTION 

The cosmological constant problem has evaded solution since its inception ||T|. If only, it 
has become more severe: the triumph of Quantum Field Theories as the correct description 
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of the fundamental interactions came at the cost of additive contributions to the cosmo- 
logical constant at big and disparate scales. It is natural to expect that the cosmological 
constant gets contributions of order Mp, where Mp is the Planck scale, from short distance 
gravitational dynamics; M^, where M\y is the mass of the W^-boson, from the phase tran- 
sition associated with electroweak symmetry breaking; Aq CD , where Aqcd is the scale of 
Quantum Chromo Dynamics, from the chiral symmetry breaking phase transition; etc. A 
proper solution to the problem has to explain how all such contributions are cancelled to 
absurdly high precision. 

An intriguing solution has been proposed f|],|3| in which spacetime is five dimensional, 
but the observable universe is constrained to a four- dimensional hypersurface, a "3-brane." 
The authors exhibit solutions to the field equations which give a flat, Poincare invariant, 
brane regardless of the value of the cosmological constant. The geometry of space includes 
naked singularities which are four dimensional hypersurfaces on either side of the brane on 
which spacetime ends. The significance of and consistency of theories with these singulari- 
ties remains unclear ||. It has been suggested || that the singularities hide a fine tuning 
equivalent to that required to set the cosmological constant to zero. In Ref. the gravita- 
tional action is modified by including a Gauss-Bonnet term, of second order in the curvature 
tensor, and it is found that the singularity can be smoothed out but only at the price of a 
fine tuning. It has also been shown Jf] that in some cases the field equations admit solutions 
that correspond to an Einstein-de Sitter universe on the brane. It has been suggested that 
this type of models may be derived from string theory ||. 

The solution is however incomplete. There is no matter in the toy model of 0,0]. It 
is necessary to incorporate matter if the solution is to be relevant to cosmology. This is 
a non-trivial issue: the standard paradigm, namely the big-bang cosmology based on a 
Friedman-Robertson- Walker (FRW) metric, is observationally very successful, so one should 
aim at reproducing, or at least approximating this paradigm once matter is added. 

We investigate inclusion of matter into these models. We discover new, exact solutions 
to the brane models with matter. They are static and therefore incompatible with observa- 
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tional cosmology. Since the solutions are not unique, it is possible that other solutions exist 
which appropriately describe the expansion of the universe. Therefore we look for time de- 
pendent solutions by linearizing the field equations about our new solutions. As will be seen, 
generically the small perturbations correspond to propagating modes or to non-expanding 
universes (static solutions). Perturbations about special, non-generic backgrounds yield 
expanding universes that seem, however, inconsistent with observational cosmology. 

The paper is organized as follows. In section |T| we briefly review the basic equations and 
the models of @|| . In section JT| we present our new solutions with non-zero matter on the 
brane. In section |V| we perform a small perturbation analysis about these new solutions. 
We present our conclusions in section |V|. 

The cosmology of brane models has been investigated in a number of papers. A general 



formulation was given in Ref. pfl. The work in Refs. [10,11] is concerned with the cosmology 



of brane models of the Randall- Sundrum type [12]. In addition, Randall- Sundrum models 



to which scalars are added have been of interest [Tj|. A method to generate solutions to 
the non-linear field equations in classes of Randall-Sundrum models with scalars was given 
in Ref. ||15|| . However, there has been little, and only very recent, work on the cosmology of 
automatically Poincare invariant branes [IB]. 



II. PRELIMINARIES 

We denote the coordinates of spacetime by x A , A = 0, ... ,4, and often use t = x° and 
y = x 4 . The 3-brane is located at y — 0. The class of spherically symmetric metrics we 
study is parameterized by three functions of t and y only 

ds 2 = G AB dx A dx B = n 2 (t, y)dt 2 - a 2 (t, y)dx 2 - b 2 (t, y)dy 2 . (2.1) 

Fixing y = we see that the metric gives a flat FRW cosmology on the brane with scale 
factor R(t') = a(t(t'), 0) where dt' = n(t, 0)dt. We will denote by g^, with v — 0, . . . , 3, 
the induced metric on the brane. 
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The models of Refs. include a scalar field 0. The action is 

S = J d 5 xVG -R+ ^(V0) 2 - Ae a * + J d 4 x [-Ve H ] . (2.2) 

We have adopted the notation of Ref. 0, save for the constants a and b which we have 
adorned with a bar to distinguish them from the metric components a 2 (t, y) and b 2 (t, y). R 
denotes the Ricci scalar. 

The peculiar normalization of the scalar field is adopted from string theory: when <p is 
a string theory dilaton its couplings are fixed. In particular, in this normalization, 6 = 2/3 
at lowest order. We are not interested solely in this particular set of string theory inspired 
parameters, so we keep the values unspecified. 

The constants A and V represent the cosmological constant in the bulk (5-dimensional 
space) and on the brane, respectively. For our analysis we set A = 0. As seen in Ref. || this 
simplifies the analysis without compromising the essential features of the model. Moreover, 
one could imagine that if the model is embedded in a supersymmetric setting, A could 
naturally vanish. The cosmological constant problem is associated with standard model 
fields' contributions to V, V ~ M% + + A^ CD H 

Einstein's equations are 

R AB _ l _ Q AB R = K 2 T AB (2 3) 

Here R AB and R are the Ricci tensor and scalar. The gravitational constant is k and 
from now on we work in units of n 2 = 1. T AB is the stress-energy tensor, which has two 
components: 

CAB 

t ab = fAB + *s(y), (2.4) 



where T AB is derived as usual by varying the action with respect to the metric, and S AB is 
a contribution from a perfect fluid of density p and pressure p on the brane, 

S A = diag (p, -p, -p, -p, 0). (2.5) 

Alternatively one may take the matter fluid to couple to the scalar, thus 



S£ = diag (e"V -e^p, -e b(p p, -e>, 0). 



(2.6) 



These two ways of writing S B correspond to distinct physical models. However the distinc- 



tion turns out to be irrelevant for the exact solutions that we present in Sec. |TJ. The field 
equation for the scalar 6 is 



-9 



5{y)bVe H = 0. 



(2.7) 



For the particular metric (2.1) Einstein's equations and the 6 field equation are 
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+ S(y) T (-Ve b ^ + p), 



(2.8) 
(2.9) 



-a 

3 



-b 2 

3 



b 2 / 
b 2 , 



%)y(^- P ), 



(2.10) 
(2.11) 



-8{y) -Vbe b 
b 



(2.12) 



Here a dot is a shorthand for d/dt and a prime for d/dy. The first four equations correspond 
to the 00, 04, 11 and 44 components of Einstein's equations. 

Conservation of the stress-energy tensor would be automatic were it derived from a 
local action. However, since a fluid component has been added on the brane, the equation 
T AB . B = contains additional non-trivial information. Conservation of energy, T 0B . B = 0, 
gives 
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while T 4B . B = yields 
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Using the field equation for the scalar, Eq. (|2.12|) , in the conservation of energy equation 
gives 

p + 3-(p + p) = (2.15) 
a 

on the brane. The second conservation equation, T 4B . B = 0, is always satisfied as a con- 
sequence of the equation of motion and Einstein's equations. Only the brane part of 
T AB . B = does not follow immediately from Eq. Q2.12| ). The identity for y = can be 
confirmed by taking the y derivative of Eq. ( |2.11| ) and then using Eqs. ( |2.8| ) and ( |2.10[ ) . 

It is stated in Ref. that for brane geometries of the form given by Eq. ( |2.1| ) the 
equation of conservation of transverse momentum, T aa -a = 0, reduces, on the brane, to 

n' a! , 
—p = 3-p. 2.16 

n a 

It is understood here that when discontinuous quantities are evaluated on the brane, like n' 
and a', they are given by their average, i.e., 

n'(y = 0) = l -{n\y = 0+) + n'(y = 0-)). (2.17) 

The equation is seldom considered any further because, for Z 2 symmetric brane-spaces, that 
is for metrics with y — > —y symmetry, which are overwhelmingly most common, the averages 
both vanish separately, a' = = n', and the equation is trivially satisfied. 

However the solutions we consider are not Z 2 symmetric. This is also true of the solutions 
in Refs. but there Eq. ( |2.16| ) is still trivially satisfied since the matter density and pressure 



both vanish. It is easy to see that our solutions below, Eqs. ( |3.8| )-( |3TT0| ), do not satisfy 
Eq. ( ^.16[ ). The reason is, in fact, that Eq. ( |2.16| ) does not apply to the case in which there are 
bulk scalars. The conservation of ^-momentum is automatically satisfied in the bulk because 
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the action is translation invariant in the bulk. So at issue here is only the conservation 
equation on the brane. Retaining only the bulk terms involving second derivatives, the 
conservation Eq. (|2.14 ) gives 



y>'$' + bS{y) 



n' a' 1 t± ( Tif a' 
—p + 3-p--Ve b n- + 3- 
n a 2 \ n a , 



0. 



(2.18) 



Following Ref. || we interpret the discontinuous derivatives on the brane as averages, and 
this will remain implicit in what follows. Integrating gives a jump equation, 



_1 A . 

36 



<V\2 



n 
n 



p -3?L p+ l V< M (VL^ 3 
a 



(2.19) 



2 \n 

We have verified that our solutions, Eqs. (|3.8j )- (|3.10| ), and also the perturbations, Eqs. (|4.7|) - 



( |1.9|) and ( f4.11| )-( [4.13j ) satisfy this equation. In fact one can prove this without reference to 
the explicit form of the solution. 



In Ref. 17 it is advocated that the correct form of the conservation of transverse mo- 



mentum equation is 



T 4 % = 0, 



(2.20) 



where the index p runs from to 3 only. However, the first term in Eq. Q2.18| ), involving 
the all important second derivative term arises from the derivative T 44 ^, which is 

ommitted from Eq. ( [2.20|) . 

Let us now describe the models studied in Refs. @||. They take p = p = 0. The solutions 
all have n(t,y) = a(t,y) and b(t,y) = 1, and are static, d = = 0. For example, case I 
studied in Ref. has A = and solutions 

;i - y/y+V + , for y > 



n = a 



(i-y/y-V~, for y < o 
¥>+iog(i -y/y+) + c, y>0 



(2.21) 
(2.22) 



<p_log(l -y/y~) + c, y<0 
The constants j± = 1/4 and \<p±\ = 3/4 are fixed by the field equations in the bulk. Case I 
has, in particular, tp± = =f3/4. The constants y± are determined by "jump" conditions, that 



is, by requiring that the second derivatives of fields in the field equations correctly reproduce 
the ^-function terms from the brane. The constant c is an irrelevant constant shift of the 
scalar field. 



III. SOLUTIONS WITH MATTER 

The models of Refs. PH provide a solution to the cosmological constant problem which 
is deficient in several ways: (1) There are naked singularities. Whether these are problematic 
remains an open question; see, for example, Refs. [fH|. (2) There is a massless scalar which 
interacts with all matter with a universal, gravity-like coupling strength. This is ruled out 
T8fl unless the coupling is made sufficiently weak. It can be arranged, however, by choosing 



the parameter b small enough. (3) It describes a static cosmology, in conflict with observation 



(see, however, Ref. |19[). (4) It does not include matter density (and pressure) on the brane. 

Here we address the last two of these problems. One hopes the two are connected: when 
matter is included in the model the universe will evolve in time. Of course, not only should 
the universe evolve, but the rate of expansion should be adequate. 

However, even after the introduction of matter the model admits static solutions. This is 
somewhat surprising, particularly if it is contrasted with the FRW cosmology which admits 
a static solution only if the matter density is precisely balanced by a cosmological term 
giving the Einstein universe. Since we are ultimately interested in time dependent solutions, 
we will look at time dependent small fluctuations about these solutions in the next section. 

We look for solutions to the field equations in the bulk with the ansatz 

a = y a , (3.1) 

n = y\ (3.2) 

b=l, (3.3) 

<f) = iflogy. (3.4) 

The field equation gives 



3a + u = l. (3.5) 

All Einstein field equations give then 

2a 2 - a + = 0. (3.6) 
9 

For definiteness we consider solutions akin to case I of Ref. 10], with 



<p = T 3 y/a/2 - a 2 , (3.7) 

where the upper and lower signs correspond to the regions y > and y < 0, respectively. 

The full solution is found by shifting y by y + {yJ) on y > 0(y < 0), and pasting these 
using the jump equations. We have 

a = A(l-^y, (3.8) 

n = N (l - , (3.9) 

<p = <p±]og(l-y/y±)+c. (3.10) 

Here A and iV are arbitrary constants that can be set to unity by a coordinate rescaling. 
The jump equations are 

Aa' b n. 



-(-Ve b * + p), (3.11) 



a 3 V2 

Aa' An' , / 1 R . \ 

2 + =6 --Fe 6 ^ + p , 3.12 

a n V 2 / 

A0' = -bbVe 1 *, (3.13) 
8 

where Aa' = a'(y = 0+) — a'(y = 0—), etc. These give three equations for five unknowns, 





at- 




11 


y+ 


y- 
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1 







Ve bc + p), (3.14) 
— - — = -p+\p + \ve bc , (3.15) 

y+ y- 3 3 



a + /2-a\ x /a^/2-a 2 _ l TTrhr 



+ -s : = -6Ke bc . (3.16) 

y+ y- 8 

These can always be solved for three unknowns (say a + and y±) in terms of p, p and two other 
unknowns (say a_ and c). The reason not all unknowns are determined is twofold. First, 
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gauge (diffeomorphism) invariance allows us to make unphysical changes to our solutions. 
This will be explained below in detail, but for now it suffices to know that one may fix the 
gauge freedom by setting, say, c = 0. And secondly, even for vanishing p and p one can find 
a class of solutions parameterized by one parameter. These are new solutions to the model 
considered in [Q, which could not be discovered there because it was assumed that n = a. 
Indeed, imposing this one has v — a which together with Eqs. ( p.5|) and ( p.7[) imply 

ot = v=\ (3.17) 

and 

<P = (3.18) 

as found in case I of 0. 

We have searched for time dependent solutions. We found one special solution, valid 
only for p = and provided 6 = 2/3. It is given by 

ds 2 = (1 - y/y ± ) 2 dt 2 - (1 - y/y ± ) 2 t 2 dx 2 - t 2 dy 2 , (3.19) 
<f> = -l]n[t(l-y/y ± )]+c, (3.20) 
— - — = -Ve 2c '\ (3.21) 

y+ y- 6 

Perplexingly, the special value 6 = 2/3 corresponds to the string theory value of this param- 
eter. Other non-static solutions were given in Ref. J7|. In this paper we will not consider 
these type of non-static matter-free solutions further. 

IV. SMALL PERTURBATIONS ABOUT LARGE MATTER DENSITY 

Armed with the new solutions with static matter density, we proceed to investigate the 
time dependence of small matter perturbations. Let us denote the static solution of the 
previous section by n , a , b , (fi , p and p . We look for solutions to the field equations, 
Eqs. (gg)-(gTg), of the form 
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n = no(l + 5n), 
a = ao(l + 5a) , 
b = b (l + 5b), 

= 00 + (50, 
p = p + $P, 

P = Po + 5p. 



(4.1) 



We count orders of the perturbative expansion parametrically in 5p and 5p. That is, we 
re-scale 5p — > e5p, count powers of e and set e = 1 at the end of the calculation. In particular 
this implies that we make no assumption as to the relative importance of temporal or spatial 



derivatives [11]. 

To derive the linearized equations in the bulk, we use again this parameterization and 
the explicit form of the zeroth order solutions, Eq. (|3.1|) . The 00, 04, 11 and 44 components 
of Einstein's equations and the field equation give 



6a" + 4-6a'--5b' + ^6<t>' = 0, 
y y 9 y 



d ( v . a r7 a r . . 4 <p , 

— —da H — 5b 5a — 5a c 

dt \y y y 9y 



0. 



25a" + -5a' + 2 on' + 5n" 

y y y 

2a + v 



-6a' 



a 



6b' + -^64>' = ^ (2<5a + 

3 y m v 

4 <p 



5n' 

y 9y 



:5a, 



1 d 



[y6<j>' + cp(35a + 5n-5b)] 



ydy - no 
The solution to these equations gives 50, 5b and 5n' in terms of 5a: 



(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 



5b 
Sn' 



-{cpSa-F), 
a 

'W-^-( 

9 a 



y 



Sa-^-( (3a-l)5a+t-F 



a |_Wq <9?/ V v 9a 
Here £ is an arbitrary constant and F is a function of t and y satisfying 



(4.7) 
(4.8) 

(4.9) 
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1 d ( dF\ 1 d 2 F _ 
y dy \ dy J n(j dt 2 

We connect the bulk solutions for y > and y < demanding continuity of the fields at 
the brane, y = 0, and using the jump equations ( |3.11| ) for the discontinuous derivatives at 
y = 0. The latter give jump conditions for the perturbations: 

A5a' = -lve Uo (5b + b5(j>) - Up - \ Po 5b, (4.11) 

A5n' = -lve Uo (5b + b6<f>) + 5p+ Up + ( Po + lp )5b, (4.12) 
o 6 6 

A5<j)' = hve l4>0 {5b + 15(f)) . (4.13) 
8 

It must be observed that both £ and F can have different values on either side of the brane 
and, moreover, that F may be discontinuous at y = 0. In addition, conservation of energy 
gives, on the brane, 

5p + 6(p +Po)5a = 5p , (4.14) 

where 5po is a constant of integration. 

There is a gauge freedom, that is, reparameterization invariance consistent with the form 
of our metric. Starting from the metric 

ds 2 = n 2 (t',y')dt' 2 - a 2 (t',y')dx 2 - b 2 (t' ,y')dy >2 } (4.15) 

we look for infinitesimal transformations 

t' = t + T(t,y) (4.16) 
y' = y + Y(t,y) (4.17) 

that leave the form of the metric invariant. Here T and Y are infinitesimal. The only 
constraint on these functions comes from the absence of off-diagonal terms in the metric: 

n 2 T'-b 2 Y = 0. (4.18) 

Under the gauge transformation the metric variations are 
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n' 

§ n = ^Y + T, (4.19) 
n 

5a = — Y" (4.20) 
a 

56 = Y\ (4.21) 

50 = <p' Y. (4.22) 

For simplicity we have indicated the variation about a static solution with bo = 1 and 
ao = ho = 0o = 0. It is instructive to check that our solutions of the field equations for 
the perturbations are invariant under these transformations, that is, that the perturbations 
Q4.19| )-( [l2~2l ) satisfy the field equations automatically. Indeed, the solution for 50, Eq. Q4.7| ), 
is satisfied provided one takes F — 0. Then, the solution for 6b, Eq. ( |4.8|) , requires £ = 0. 
Finally the solution for 5n', Eq. Q4.9|) , is satisfied provided 

T' — \y — 0, 

which is a consequence of the condition ( 4.18| ). 



One may fix the gauge by imposing, for example, 

6b{y, t) = and 5<p(y = 0, t) = 0. (4.23) 

There is some residual gauge freedom: there are further transformations with Y = and 
T = T(t). These are uninteresting time reparameterizations. Were we to impose a gauge 
condition on any of the fields that have a discontinuous derivative on the brane, we could only 
require that it vansihes for either y > or y < 0, but not both, since the reparametrization 
functions are smooth. 

We are ready to present our solution to the linearized field equations. The solutions to 
the bulk equations express 50, 5b and 5n' in terms of 5a and F. Once F is determined in 
terms of 5a our task is to determine 5a only. We will proceed as follows. First we use the 
gauge conditions and our solutions in the bulk, Eqs. (|4.7|)-(fO|), to eliminate the function 
F everywhere, and to relate 5a' to 5a on the brane. Then we consider the jump equations 
which impose further restrictions on 5a and its derivatives on the brane. These constraints 
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on the brane are used, finally, to determine 8a in the bulk. The cosmlogy depends only on 
the fields on the brane. Therefore, we concentrate on determining as fully as possible the 
fields on the brane. 

We eliminate the function F by fixing the gauge as in Eq. ( |4.23| ) and using our solution 
of the field equations in the bulk, Eq. Q4.8p, thus 



F = -- [(ySa)' - . (4.24) 

4 ip 

We can now fix 8a! on the brane in terms of 8a. Consider the constraints from requiring 
continuity of 5<fi at y = 0. We have in addition the gauge choice in Eq. ( f4.23| ), 5<p(y — 0, t) — 
0, so we obtain 

(4a + -l)5a = (y + 8a' + + Z + ), (4.25) 
(4a_ - l)8a = (y^5a'_ + £_). (4.26) 

Here and below we denote the limiting values of fields as y — > 0± with a corresponding 
subscript, eg, 5a' ± = \im y _ ±5a' . 

The first jump equation, Eq. ( |4.11|) , gives a constraint between the constants of integra- 



tion. From the jump equations for the exact solution, Eqs. ( |3.14|) and Q3.15|) , one has 



4a,-l 4a_-l . A . 
Po+Po = — + • 4.27 

y+ y- 



Using this and Eqs. ([4.25|) and ( |4.26| ) in Eq. ( |4.11|) we obtain 



^Pb = — + — • (4-28) 

3 y+ y- 



The remaining two jump equations, ( |4.12j ) and ( |4.13| ), involve 5a± in addition to 5a' ± , 5a 



and 5a. They can be solved to give, on the brane, 5a± in terms of 5a and 5a. However, 
the resulting expresions are long and not terribly illuminating so we refrain form presenting 
them here. 

The jump equations do not fix the time behavior of 5a on the brane. We now describe 
a procedure that determines the time behavior and the bulk dependence of 5a. Note that 
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since F is given in terms of 5a through Eq. (|4.24j ), one must now impose that 5a satisfy the 
wave-like equation ( |4.10| ). The solution to Eq. ( |4.10| ) is straightforward, 



F ± (y,t) = / d»**F ± M IT^JL, ■ (4-29) 

^ 3a± I I / 

where Jo is a Bessel function. There is a second solution involving the Neuman function. 
However, we have dismissed it since it is arbitrarily large as \y— y±\ — > 0, which is outside the 
validity of perturbation theory. The Fourier coefficient F±(u) is determined by the function 
on the brane, F(0, t). Using the conditions on derivatives of 5a on the brane one may express 
this solely in terms of 5a: 

9a±(l - 2a±) 



-6a{0,t) = J due iujt F±{u). (4.30) 



2 ip ± 

Thus, knowledge of 5a on the brane determines all fields everywhere. The time- 
dependence of 5a, however, is not arbitrary. Consistency of Eq. ( |4.29j ) and the jump equa- 
tions gives an integro-differential equation for 5a on the brane. To see this, differentiate 
Eq. ( f4.29| ) and evaluate at y — 0: 

J ( — \v \^ a± ) 

- y ± 5al + 25a' ± = [ due^(4a ± - 2)5~a(u)u\y ± \^- 1 1 ^ ± (4.31) 

where (1 — 2a±)5a{u)) = 2<p±/9a>±F±(uj). Now the left hand side of this equation can be 
expressed in terms of 5a and 5a through the jump equations. One obtains 

- y±5al + 25a' ± = A±5a + B±5a + C±, (4.32) 

where A± and B± are constants given in terms of the parameters of the background solution, 
and C± also depends on these parameters but is, in addition, linear in the small constants 
£ ± . We have ommitted the complicated and long expressions for these. Thus, we obtain an 
equation that determines the time dependence of 5a 



A particular solution to this eqution is easily found by Fourier transform 



A±5a + B±5a + C± = J du e iuJt {4a± - 2)5a{u)uj\y±\ 3a± - 1 j 3 ^ - (4.33) 
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Sa(Q,t) 



C±/B ± . 



(4.34) 



To this solution one may add arbitrary linear combinations of solutions to the associated 
homogeneous equation (obtained by setting C± = 0). These solutions have time dependence 
of the form of simple exponentials, expiuj Q t, with characterisitic frequencies u that are 
solutions to 



There is an infinite number of solutions to this equation. The function Jo has an infinite 
number of simple zeros. Between two succesive zeros the last term takes on any value. There 
may be additional complex solutions when A±B± < 0. 

The physical interpretation of our solution is straightforward. The infinite set of oscila- 
tory modes simply correspond to field excitations in the bulk. Neither their amplitudes nor 
their frequencies depend on the additional density perturbation 5po- This is to be expected 
because even in the absence of additional matter there can be propagating gravitational and 
scalar field waves. It is only the particular solution that actually depends on the added 
matter. Ignoring the possible excitation of propagating modes, we see that the new solution 
is static. 

There are some caveats. This conclusion does not hold if the parameters of the back- 
ground are such that B± = or if there are complex solutions ujq. In the former case, if in 
addition C± = 0, there are further solutions of the form 



This is of particular interest beacuse the perturbations may now be time dependent. In 
fact, such time dependent perturbations are expected since they must correspond to the 



where A is an arbitrary constant. We have verified that B± = and A± ^ for parameters 
corresponding to the background of Ref. J7| (a± = 1/4, ip± = =F3/4, p Q = 5po = 0). There 
are also solutions of the form 5a cx t 2 . If there are modes of complex Uq, the solutions are not 




(4.35) 



5a oc t. 



(4.36) 



linearization of the exact solutions of Ref. |7J for which a oc exp(/(?/)+V / At) — l+f(y) + VXt, 
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simply oscilatory, but will in addition have exponential dependence. None of these solutions 
are suitable for observational cosmology. 

V. CONCLUSIONS 

The proposed solution to the cosmological constant problem of Refs. |2|||] is incomplete in 
that it does not include matter on the brane, i.e., in our universe. We have found new exact 
solutions including matter on the brane. These solutions are static and therefore describe 
an unacceptable cosmology. 

However, the equations admit other solutions even under the same assumptions on the 
symmetry of the metric. This had been recognized in Ref. which found curved brane 
solutions to the model dubbed case I in Ref. [||. Here we have studied new solutions 
obtained as small perturbations about our new static solutions with matter. The small 
perturbations that are not simply propagating waves are generally static. In special cases 
the non-propagating, small perturbations grow linearly with time, but we have identified 
these as the De-Sitter-like solutions of Ref. ]/]]. On this basis it is tempting to rule out these 
as viable cosmologies. 

In standard FRW cosmology the evolution of the scale factor is completely determined 
once the equation of state is fixed. Once spherical symmetry is chosen and matter is specified, 
Einstein's equations determine the cosmology. However, this is not the case in the peculiar 
brane cosmologies of Refs. 0|§. To understand what is happening one could consider this 
as an initial value problem. In standard cosmology if the metric and matter content were 
specified at an initial time (in a fixed gauge), one could evolve forward using the field 
equations. Thus one would recover the standard picture. Clearly this is not the case of the 
brane models. What else must be specified and why? This is an interesting question that 
we hope to explore further. Our guess is that the naked singularities introduce additional 
information that has been implicitly specified. In the absence of a new general principle that 
specifies these additional data one would have to give up the notion of causality (at least 
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on a global scale). This may be the price one must pay in order to solve the cosmological 
constant problem. 
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